Modeling Copper Cables and Connectors* 
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High-speed data busses constructed from copper components are both difficult to model and 
important to accurately simulate. Using analytic signal theory and an extension of the skin effect 
for non-uniform conductivity, experimental data taken from IEEE f 394- f 995 cables and connectors is 
used to predict the transfer function of a 4.5 meter cable with two connectors. With these techniques 
any length of cable between two connectors can be simulated with high accuracy, including all 
the cable and connector losses and the infinitely summed standing waves that occur between the 
connectors. 

Keywords: copper interconnect, analytic signal, skin effect, transmission line 



I. INTRODUCTION 

Modern high-speed busses send data streams across 
long runs of copper, often at data rates in excess of 1 
gigabit/second. Connectors at the ends of copper cables 
create discontinuities that phase shift the electrical sig- 
nals and cause reflections. Mathematically characterizing 
these copper systems with high accuracy is challenging 
and often performed with awkward lumped-circuit mod- 
els. In this paper I demonstrate methods that provide a 
new level of accuracy and convenience in modeling such 
systems. To demonstrate the accuracy of this method the 
IEEE 1394-1995 connector and cable system is analyzed 
for jitter and losses. 



be treated as a transmission line with a simple skin depth 
effect - real wires have non-uniform conductivity due to 
stresses and alloying effects in their manufacturing. Fi- 
nally, the entire system acts as a lossy resonator with 
standing wave modes. 

Accurate modeling and simulation of this simple a sys- 
tem is both difficult and necessary. High-speed cable sys- 
tems may have frequency-dependent losses specified to 
the 1 dB level and jitter budgets defined to 100 ps - for 
5 meters of cable. Errors in the predicted waveforms are 
catastrophic in systems that may ship in the 100 million 
units per year volume. 



III. ANALYTIC SIGNALS 
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II. COPPER INTERCONNECT 

Figure |l] shows a simple copper interconnect consisting 
of two connectors and a long run of two parallel copper 
wires. 

We restrict the discussion to the differential mode of 
excitation, excluding common mode effects - and also 
neglect cross-talk and EMI issues. The connectors have 
complex effects upon the phase and amplitude of each 
Fourier component of the electrical waves incident upon 
them. They cannot be treated as simple lumped ele- 
ments in high-speed systems. Next, the wires have dis- 
tributed effects upon the electrical waves and should not 



FIG. 1: Copper interconnect 



The analytic signal Q is a complex- number represen- 
tation of a signal. The basic notion is that we start with 
a real- valued signal which has no explicit phase informa- 
tion, make this into the corresponding analytic signal, 
and thereby determine the phase factors of the complex 
signal. (See Appendix [X]). Analytic signals are essential 
for proper modeling of copper interconnects. 



A. Amplitude-Phase-Frequency ambiguity 



Following Q we write a signal in the form 



u(t) — a(t) cos <f>(t) 
= a(t) cos(ui(t)t 
u(t) € K 



•*(*)), 



(1) 
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a(t), and w(t) = d<p/dt are the amplitude, phase, 

and frequency of the signal. However, for a real-valued 
signal u(t), this is an equation relating two unknowns to 
a single known. It may be argued that given the smooth 
envelope of maximums and the zero-crossings we know 
what the amplitude and frequency are. To a certain ex- 
tent this viewpoint is effective for narrowband signals. It 
is totally inapplicable for wideband signals. 
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It is advantageous to work with complex signals 

w(t) = u{t) + iv{t) = a(i)e^ (t) (2) 

constructed by adding an imaginary part v(t) to the real 
signal u(t). For complex signals the amplitude, phase 
and fequency are well defined 



a 2 (t) = u 2 + v 2 = \wf 



4>{t) = arctan[— ] 



(3) 



(4) 



The continuous-time analytic signal z(t) corresponding 
to x(t) is most simply defined in the frequency domain 
as 




(8) 



which is inverse transformed to obtain z(t). Note that the 
value of Z(f) at / = is defined to reflect the mathemat- 
ical behavior of the Fourier transform at a discontinuity, 
which yields the average of the values on either side of the 
discontinuity. It can be shown that the imaginary part 
of z(t) is the Hilbert transform of the real signal x(t). 



uj(t) = Im[— ] 
w 



(5) 



If v(t) is the Hilbert transform of u(t) then w(t) is the 
analytic signal. The Hilbert transform is 



H[u(t)} = 



1 



t - s~ 



ds 



(6) 



VakmanQ has demonstrated that only the Hilbcrt- 
transform imaginary part of the complex signal satisfies 
the following conditions: 1) amplitude continuity and dif- 
ferentiability under small perturbations, 2) phase-value 
independence of signal scaling, and, 3) harmonic corre- 
spondence, i.e., sine waves of constant frequency and am- 
plitude retain their values. 

The analytic signal provides a cogent and global defi- 
nition of the amplitude, phase and frequency of a signal. 
The analytic signal is implicitly implemented in radio 
devices such as filters and modulators. Finally, the ana- 
lytic signal provides accurate estimates of the amplitude, 
phase and frequency of wide-band signals even in the 
presence of substantial noise. 



B. Fourier transform for analytic signals 

To efficiently convert time-dependent real signals to 
analytic signals we use the Fourier transform. 



1. Continuous-time analytic signal 

Following Marple|| let x(t) be a real-valued, finite- 
energy signal defined over the temporal interval — oo < 
t < oo with the corresponding continuous-time Fourier 
transform (CTFT) 



*(/) 



x(t)e^ i2jTft Ut 



(J) 



defined over the frequency interval — oo < / < oo. Be- 
cause x(t) is real, the CFTF is complex-conjugate sym- 
metric, X(-f)=X*(f). 



2. Discrete-Time "Analytic" signal properties 

Again following Marplejj| we consider the properties 
appropriate for an analytic- like discrete-time signal z[n] 
corresponding to a real- valued discrete-time signal x[n] 
of finite duration T. 

The spectrum of the signal x[n] is obtained from the 
discrete-time Fourier transform (DTFT) 



N-l 

X(f) = T ^ x[n]e^ l27rfnT \ 
n=0 

/ = 0,l,...,iV-l 



(9) 



The DTFT is evaluated with a fast Fourier transform. 
Next form the TV-point one-sided discrete-time " analytic" 
signal transform 



Z(m) 



X(Q), 
2X(m), 



m 





N 



1 < TO < f - 1 



(10) 



0, f + l<m<N-l 
Compute using an TV-point inverse DTFT 



z\n\ 



N-l 

— y zMe (+427rmn/JV) 

m=0 



(11) 



to yield the complex discrete-time "analytic" signal for 
the same sample rate as the original signal x[n]. 

One justification for setting the special value at m — 
Y is that Hz and the Nyquist frequency terms are 
shared boundaries between negative and positive fre- 
quency halves of the periodic spectrum, and division 
into respective one-sided positive and negative spectra 
requires that these terms be split between the two spec- 
tral halves. 



IV. EXPERIMENTAL METHODS 

We set up in the laboratory a standard IEEE 1394-1995 
connector and cable assembly which was interfaced to a 
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Tektronics Digital Sampling Scope with a Time Domain 
Reflectometry module (TDR). TDR data from produc- 
tion cables and connectors was analyzed for impedance 
discontinuities and then used to model data-dependent 
jitter. 



A. Time domain renectometry 

TDR directly measures the reflected voltage from a 
test fixture as a function of delay time. A step function 
with a very short risetime is launched into the connected 
system. After the launch trigger the input of the system 
is monitored for voltage levels. Voltage levels are affected 
by spatial variations in the impedance of the connected 
system, as changes in impedance scatter the outgoing 
waveform back towards the source point. 



B. Digital sampling oscilloscope operation 

The Tektronix 11801B oscilloscope contains a two- 
channel 20GHz bandwidth sampling head. The oscillo- 
scope functions in equivalent-time sampling. A sampling 
interval is defined as the horizontal time divided by the 
record length (number of defined samples) . As an exam- 
ple, for a sample interval of 200ps/div the oscilloscope 
advances by one sample interval (200ps) after each trig- 
ger event. Hence, for a record length of 512, it would 
take 512 trigger events for the oscilloscope to advance to 
the end of the screen. 

TDR operation adds a step recovery diode circuit that 
generates a step pulse with a 40 ps rise time into the 
Device Under Test at the time strobe. All the reflected 
signals are sampled by the DSO system to show the time 
profile of the DUT's reflections. 



V. EXPERIMENTAL WORK 

The connector was a surface-mount style and the cable 
was a D7 style. Only one connector and cable assembly 
was tested though our results correspond closely to data 
given to us confidentially by other companies. 

The Figure |^ shows the differential TDR signal for 
a production 1394-1995 connector and cable. Figure || 
shows the TDR data for a connector with the cable re- 
moved. We see there is a substantial dip in the connec- 
tor's impedance. 

It is useful to convert TDR data into an equivalent 
impedance profile. For differential TDR the conversion 
formula is: 



P.ef lectometry profile 



Z = 100 ft ( 



1 



(12) 



Note that the cable impedance in Figure [| is above the 
100 Ohms level (about 105 Ohms). The ideal 1394 cable 
differential impedance is 110 Ohms. 
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FIG. 2: TDR of connector and cable 
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FIG. 3: TDR of connector with cable disconnected 



A. Complex impedance of the connector 

Following the numerical algorithm given above we now 
create the analytic signal representation of Figure ||'s 
impedance profile (which has no phase information). 
This algorithm will generate the imaginary component of 
the impedance. Figure ^| is the resulting analytic signal. 
The top line is the real part of the analytic signal and the 
bottom, broken line, is the imaginary part. Taken as a 
complex function this diagram is the complex impedance 
profile of the connector. A lumped-circuit approximation 
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FIG. 4: Equivalent impedance profile 
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Analytic Signal o± Impedance 
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FIG. 6: Synthetic pulse 



FIG. 5: Equivalent impedance profile 



for this would require a series of capacitors and inductors 
to match each minimum and maximum in the trace - a 
large number of components. Additionally, such a se- 
ries of components would have secondary reflections that 
would need to be added into the fitting procedure so that 
the total, consistent, electrical system matched the TDR 
data. 



B. S-parameters 

We now want to derive the frequency-dependent S- 
parameters for transmission through and reflectance off 
the connector. 

The definition of the transmission coefficient, S12, is: 



Reflected Pulse 



a u A _ TV, ,\ _ 2Z loadM 

Z load( w ) + Z o(w) 



(13) 



Similarly the definition of the reflection coefficient Sn is 

<7 (,,\- R(< A - Z load(^)" Z oH 
z loadM + Mu) 



(14) 



The most straightforward way to derive these coeffi- 
cients for the connector is to create a synthetic pulse with 
an exactly known Fourier transform (see Figure [3]), use 
the time profile of the TDR derived connector complex 
impedance to calculate reflectance (see Figure |7] - mag- 
nitude only) and transmission (see Figure M - magnitude 
only) synthetic pulse patterns, then Fourier transform 
the time profile of these pulses to get the analytic spec- 
tra. One subtlety is to use the variation of the measured 
impedance from the nominal 11 Oil value of the trans- 
mission line, not the absolute impedance. It is only the 
variations from the uniform value that scatters waves into 
the backward direction. For accuracy we use the actual 
transmission line impedance value measured (by suitable 
averaging away from the connector). A further subtlety 
is that we have implicitly summed up all the internal re- 
flections inside the connector by using the TDR data. 



0.5 1. \ v /1.51M 2 V 2.5 3. 





FIG. 7: Reflected pulse 



Finally, we divide the analytic spectrum of the re- 
flected and transmitted pulses by the analytic spectrum 
of the original pulse to get the complex ratio at each 
frequency, which is the transfer ratio as a function of fre- 
quency. This procedure simulates using a Network Ana- 
Iyer. Figure |9| is the S matrix function for S12, Figure 
|l0| is Sn. Using the above complex- valued reflection 
function we can validate this approach by simulating the 
reflection of an ideal step function. This should match 
the original p data. Figure [ll] shows an almost exact fit 
with the experimental data (dashed and offset by 0.05 
from the simulation). 



VI. CABLES 

The next element to consider for copper interconnect 
is the cable. We need to determine an accurate model 
that captures the experimental data. 
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FIG. 8: Transmitted pulse 
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FIG. 12: TDR data for 4.5 meter cable 
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outward trip and on the return trip. By analyzing the 
amplitude and phase changes long cables apply to each 
Fourier component we can derive a very accurate model 
of the cable's distributed parameters, e.g., loss per unit 
length, incremental inductance, etc. 



B. Cable equations 



FIG. 10: | Sul 

A. Experimental data 

Figure [l2| shows the experimental data for a 4.5 me- 
ter D7 cable with a shorted end. This data was taken 
by using TDR. The time axis measures roundtrip delay 
times. The ringing waveform on the left side is due to the 
fixture used to couple the TDR equipment to the cable 
and is excised from successive analysis. There are several 
things to note about this figure. First there is a linear 
rise in resistance as time progresses. This is due to the 
DC resistance of the cable. The TDR pulse scatters off 
this incremental resistance and produces a gradual series 
of return signals that add up to the observed ramp. The 
next feature is the abrupt drop at about 50 ns time delay. 
This is where the TDR pulse hits the short at the end of 
the cable. For a very short cable this drop would be al- 
most perfectly vertical. For long cables the return pulse 
has been modified by the cable imperfections both on the 
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Two-conductor cables operating purely in differential 
mode are described by the partial differential equations 

d z v(z,t) — —ri(z,t) — ld t i(z,t) (15) 



d z i(z, t) = — gv(z, t) — cdtv(z, t) (16) 

where r is the resistance of the cable per unit length, c is 
the capacitance, and / is the inductance. More generally 
we can make r frequency dependent, g is the transverse 
conductance between the two conductors, the g term is 
neglected by assuming there are no dielectric losses. As 
a subtle point, c and / are assumed to be determined by 
the cable's wire diameters, i.e., we are assuming perfect 
conductors and perfect shielding in the above partial dif- 
ferential equations. We will correct for the effects of skin 
depth on these equations later. 

We need to solve these equations for one-way energy 
flow on an infinitely long cable. First we assume a har- 
monic time dependence 

v'{z)e~ ttu = - (ri{z)e- ltuJ ) +ilcji(z)e- ltuJ (17) 



i\z)e- lt " =icujv{z)e- ltw (18) 

Removing the common time dependence 

v [z) = — (r i(z)) + tl w i(z) , 

i'(z) = icluv(z) (19) 



FIG. 11: Simulated TDR vs real TDR 



Now we solve these differential equations 
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f l + e 2zV-icru>-clu>A f 1 + e 2zV-tcr^-cl^\ ( r _ ^ w ) C(2) 

"(*) = > ; — (2°) 

2 e zV-icru-cI^ JgzV-icru-cl^ y/-tcrUJ - cliO 2 



5 c (-l + e 22V -" : ™- d ' 1,2 j wC(l) (l + e 2zV - ,cr ^ clu! j C(2) 



i(-s) = , ., ' + - j (21) 

g z v — i crw — cl u) i c V LO cl UJ 2 2 C z ^ ~ 1 c r u ~ c ^ 



Collecting coefBcients between outgoing and ingoing waves 



-^cuC(\) C(2) , 

gW-icr^-ci^ ' VV-Jcrw - cIlo 2 2 J 

C(l) + rC(2) _ j^C(2) 
V(Z) = === (23) 

g3 y/ — l CT UJ — cl UJ 

+ e W -»c^- c ^ ( C(l) rC(2) + \luC{2) 



We now define C(2) to eliminate the incoming wave 



C(2) 



y—icruj — cluj 2 



We have 



;(«) = e ^V-( c "(»'+^))c(l) 



zee 



z V cr uj — cl l 



»C(1) 



\/ —icruj — clt 



(24) 



(25) 



(26) 



We now derive the cable's attenuation factor. Focusing 
upon v(z) and setting C(l) = 1 we have 

v(z) = e z \M C "("+M) (27) 

To further simplify this equation we note that r is a small 
impedance compared to the inductive impedance of the 
cable, for all high enough frequencies. By expanding the 
square root around the large inductive term and dropping 
the imaginary phase-determining part, we derive the loss 
per unit length of the cable. 



loss — e 2 vT 



In log terms 



logLoss 



- (Vcr z) 
2~Tl 



Which results in 



-2 Vl logLoss 



(28) 



(29) 



(30) 



This formula allows us to determine r as a function of 
frequency from the experimentally measured cable loss 
versus frequency relationship. We set z to the measured 
length of the cable, times 2 to account for the roundtrip, 
to derive the r per unit length values. This frequency- 
dependent r function produces a much more accurate 
model for real cables then a simple skin depth effect 
model. 

In Appendix |b] we show that the skin effect produces 
both a resistive term, due to losses in the wire from the 
finite conductivity of the metal, and an additional induc- 
tive term, due to the finite penetration of magnetic flux 
into the wire. After we measure the real part of r from 
the experimental cable data we add back the skin effect's 
additional inductive effects by using the " complex" resis- 
tance value: r extracted (l - i). 



C. Experimental cable results 

Figure |lj shows the full experimental data set. The 
most interesting part of the data set is the sharp return 
step. To isolate just this element we window the data 
set with a gaussian envelope, centered upon the step, 
but broad enough to preserve the lower frequency com- 
ponents. Figure [13] shows the result. 

To compare the experimental data with a known wave- 
form we do the same gaussian windowing with an ideal 
return pulse that has been scaled and positioned to best 
match the experimental data. Figure [l4| shows this wave- 
form. 

After Fourier transforming both of these windowed 
pulses we divide the experimental Fourier components 
by the ideal Fourier components, frequency by frequency. 
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FIG. 13: Windowed cable data 



FIG. 16: Cable resistance 
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FIG. 14: Windowed ideal model 




FIG. 17: Experimental edge (bold) compared with the mod- 
eled edge 



Figure |lj shows the loss function of the measured cable. 

Using the formula above we derive the cable's mea- 
sured resistance factor shown in Figure [l6|. 

Figure [l?] shows a close-up of the experimental TDR's 
step edge versus a simulation based upon the mathemat- 
ical methods we've developed above. The bold plot is the 
data. We see reasonable agreement. 



H(lo) is the cable transfer function, T(uS) is the trans- 
mission coefficient for the connectors, and R(uo) is the 
reflection coefficient for the connectors. 

Using the connector and cable models we've derived we 
have the composite transfer funcition shown in Figure 
|lS| , where we've made the cable very short. Signficant 
resonances are visible. 



VII. CONNECTOR PLUS CABLE SYSTEMS 



VIII. EYE DIAGRAMS 



Given two connectors and a length of connecting cable 
we have an infinity of reflections which sum up to a net 
transfer function defined by 



H{lo)T 2 {lo) 
1-R 2 (uj)H 2 (l 



(31) 



An important quantity to determine for a copper in- 
terconnect system is the eye diagram. We synthesize a 
stressful set of data transitions, apply them to the copper 
system transfer function, and see how the output pattern 
looks. The basic technique is to convert a bit pattern to 
a pulse pattern with the right timing, e.g., rise/fall times, 
etc. We produce a discrete, time-sampled waveform from 



1394 Cable Differential Mode Losses 
(4.5 meters ) 
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FIG. 15: Cable losses 



FIG. 18: Two connectors and a short length of cable 
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400 Hbs @ 4.5 meters 
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FIG. 19: 400 Mbps eye diagram 



the pulse pattern, Fourier transform it, multiply by the 
transfer function, then inverse Fourier transform to get 
the output waveform in time. 

We use the following bit pattern to test the ca- 
ble/connectors system: 11000001 01001111 10101000 
00000000 00000000 00000000 01011111 11111111 
11111111 11111110 • • • repeating. 

Figure [l^ shows the result of this analysis. We have 
sliced the output waveform into 1 bit-time length seg- 
ments and then aligned and superimposed these slices. 
The result of the simulations is a very accurate model 
of experimentally determined eye diagrams (not shown 
here) . Eye diagram simulations permit the determination 
of useful things like total system margins, jitter budgets, 
etc. 



IX. CONCLUSIONS 

TDR Oscilloscopes are readily available pieces of 
equipment. By using TDR techniques along with the 
mathematical methods described in this paper it is pos- 
sible to characterize and model many of the interesting 
features of a copper interconnect system. 



APPENDIX A: THE ANALYTIC SIGNAL 

We have used analytic signals (see GaborjlJ) to de- 
termine the phases of complex- valued transfer functions 
when all we originally had was the real part of the func- 
tion. In this appendix the theory behind analytic signals 
is discussed. 

The essential assumption behind analytic signals is 
that we are dealing with a linear system that obeys 
causality, i.e., responses occur after inputs in time. The 
most general way to express this is by the integral equa- 
tion: 



dtK(t-t)g(t) 



(Al) 



K(t—t) is the Green's function for the system (e.g., cable, 
connector), g(t) in the input (e.g., incident voltage wave- 
form), and 4>(t) is the system's response (e.g., reflected, 




FIG. 20: Integration contour 



transmitted waves). This integral equation comprehends 
all the time lags of the linear system. We require causal- 
ity which means that K(t) = for t < 0. The function of 
interest for our system is the Fourier transform of K(t) 
and the condition Kit) = for t < implies that the 
Fourier transform 



f(z) = / dte Uz K(t) 



(A2) 



of K(t) has no singularities for Im(z) > 0, and that 
f(z) —>■ for \z\ — ► oo in the upper half z-plane. 

We now apply Cauchy's formula to a contour consisting 
of the real axis and a large upper semicircle, see Figure 
|20| . By the above condition on f(z) we have the con- 
tribution from the semicircle to the line integral along 
the contour going to zero. This implies that Cauchy's 
formula reduces to 



/(*) = 5?- / 



/(C) 



real axis C 



d^, Imz > 



(A3) 



The function of physical interest is f{z) on the real axis 
(i.e., real- valued frequencies). Define the function F(x) 
to be 



F (x) = lim f(x + ie) 

e— >0 

where e is positive. Then it follows that 

F(x) 



(A4) 



2inF(x) 



lim 



= P 



F(x) 



dx 



x — x 



dx + inF(x) 



(A5) 



where P designates the Cauchy principal value. Thus 

F(x) = —P [ -P^-df (A6) 

«7T J_ t 



1 f C 

ReF(x) = -P / 

7T J-> 



X — X 



ImF(f) 



x — x 



dx, 



(A7) 
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ImF(x) = — P 

IT 



ReF(x) 



dx, 



(A8) 



These are dispersion relations for the function of phys- 
ical interest F(x) and the last equation clearly shows 
how to calculate the imaginary component of any trans- 
fer function from knowledge of its real part. This equa- 
tion is the Hilbert transform. Note in our application x 
is frequency. 

As a simple example let us consider a simple low-pass 
RC circuit with unit R and unit C in series. Vj n is ap- 
plied across the series combination and V ou i is measured 
across the capacitor. The complex transfer function is 



1 



V in 1 - J- 



(A9) 



which can be expressed in real and imaginary parts ex- 
plicitly as 



V m 



1 



Vt, 



l+LO 2 1+LU 2 



real part 



1 



(A10) 



(All) 



We now assume that all we know is the real part of 
the transfer function. The goal is to use analytic signal 
theory to derive the correct imaginary part. Negative 
frequencies don't make sense physically in this situation 
but this form of F possesses the symmetry: 



/(-*) = +/* CO 



where * denotes complex conjugation. This implies the 
relationship 



i r F 

Fix) = —P - 

ITT X 



1 f° F*(x) 
D ' y 'dx 



X + X 



1 

— p 

VK 



For the imaginary part 
ImF(x) 



^P 



oo 



F(x) 
ax 



-ReF(x) 

-dx 



[I X "T X 

00 -ReF(x) , , 

; dX 

Q X X 



2_ p r***m d . 



7T ./,. ./•- — X 2 



which simplifies to 
ImF(x) 

Putting in the explicit value of F 

-2 f°° 
ImF(w) = — P / - 

7T Jo U 



rdu) 



-2 



-P 



(A13) 



(A14) 



(A15) 



7T J Q (Lb 2 - U0 2 ){1 +UJ 2 ) 



dib 



(A16) 



(A12) Given the definition of the Principal part: 



P 



(LU 2 -LU 2 )(l+LU 2 ) 



deb 



lim 

5^0 



(LU 2 - U 2 ){1 + Lb 2 ) 

I 



dub 



UJ + S 



(W 2 - UJ 2 )(l + Lb 2 ) 



dib 



(A17) 



where 5 is positive. Performing the integrations 
Im_F(a;) 



1 



(A18) 



which is exactly the imaginary part of the transfer func- 
tion. 



APPENDIX B: SKIN EFFECT 

Accurate prediction of the behavior of copper trans- 
mission lines requires proper modeling of the losses and 
incremental inductance factors caused by the skin effect. 
Here we treat the skin effect without the usual assump- 
tion of spatially constant conductivity. 

Following JacksonQwe restrict Maxwell's equations to 
weak fields, linear, isotropic-response media, and the 



quasi-static regime - where we neglect the displacement 
current contribution to the magnetic field. 

V X H = J 
V • B = 

V X E + = 



J = aE 

B = /xH 



(Bl) 



With the definition of the vector potential, 
B — V X A, Faraday's law shows that the curl of 
E + dA/dt vanishes. This implies that we can write 
E — —dA/dt — V$. Assuming negligible free charge 
and that the time-varying B is the sole source of the 
electric field, we may set the scalar potential $ to 
zero and have E — —dA/dt. Note that we have the 
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subsidiary conditions V • E — and V • A — 0. For 

media of uniform, frequency-independent permeability 
ji, Ampere's law can be written V X B —pj — paE. 
Elimination of B and E in favor of A and use of the 
vector identity, V X V X A — V(V • A) - V 2 A, 
yields the diffusion equation for the vector potential 

V 2 A = pa°— (B2) 

This equation, which obviously also holds for the elec- 
tric field is valid for spatially varying, but frequency- 
independent conductivity. If the conductivity is constant 
in space it follows that the magnetic induction and the 
current density also satisfy the same diffusion equation. 
This situation gives the usual skin depth effect. Constant 
conductivity is a poor approximation for real wires which 
are subject to varying alloy concentrations and stresses 
due to processing. The goal of this appendix is to de- 
rive the overall relationship between the driving current 
and electric fields in realistic wires with spatially varying 
conductivity. 

To proceed we assume cylindrical coordinates along the 
wire, harmonic time dependence, and use complex phasor 
notation, A(p, <f>, z,t) — A(p, <fi, z)e~ luJt . After dropping 
the common exponential factor 

V 2 A(p, <t>, z) = -iupa(p, <f), z)A(p, </>, z) (B3) 



Assuming the wire diameter is much smaller that the 
wavelengths of interest, we drop the <fr dependence. Given 
that a constant amount of current is flowing along the z 
axis of the wire we see that the curl of A at the surface of 
the wire must be constant, independent of z. This follows 
from the relationship between B and A and Gauss's law 
relating the magnetic field around a contour and the en- 
closed current. E is very small in the deep interior of the 
wire at high enough frequencies, thus A is very small at 
the center of the wire. We also assume that E is parallel 
to the wire at its surface, given there is no free charge on 
the surface. The only non-trivial dependence is that of 
a and A on p and z. We ignore ct's and j4's dependence 
upon z since each different position along the wire will 
have an independent voltage drop due to the very small 
skin depth expected at interesting frequencies - the net 
voltage to driving current ratio for the entire cable is 
simply the average along the wire's length. Thus we as- 
sume a uniform and constant dependence along z for the 
diffusion equation. We now have a set of equations to 
solve 



a , \ I d , dA z , ,d 2 A z 
V 2 A z (p,z) = - — (p—^) + (- 



p dp dp 



dz 2 



0) = -iajpa(p)A z (p), 



(B4) 



V 2 A p {p,z)= l -^-{p d 4^) + ( d2Ap 



p dp dp 



dz 2 



0) = -iujpa{p)A p {p), 



(B5) 



o)-^M)i 

U > dp ' lp=r wire" 



B 



surface 



r 



(B6) 



4,(0) = 4,(0) so, 



(B7) 



Mr w ive) = (B8) 

4( r wire) = ~ since the electric field at the surface of 
the wire is strictly oriented along the wire's axis. Given 
the boundary conditions we see the A p always remains 
zero. There are no terms in the partial differential equa- 
tion for A p that can drive a non-zero value anywhere, A p 
being decoupled from A z . All that remains is 

1 d dA 

V 2 A z (p, z) = -_(p_£) = -ioopa(p)A z (p), (B9) 



p dp dp 

dp wire 



B 



surface > 



(B10) 



4(0) = o 



(Bll) 



As a first step we solve these equations for spatially 
constant a 



A z (p) = Jo(- \h-Y^/p,oup)c\ + Y (-V-lv / pcrLup)c2 

(B12) 

Y (0) = — oo, thus to preserve the small value of A z at 
the origen of the wire we set C2 to zero. 



dA z (p) 

dp 



(B13) 



= Ji (- y^iy/ji^jp)^ y^iy/ji^; (Bi4) 
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Setting the derivative at p = r : „ equal to the surface magnetic field we have the full solution 



Mp) = - 



B. 



surface 



l^fpaujJi (- ^l^/JIaUJr wive ) 



(B15) 



dA z (fi) 

dp 



-^surface 



Ji(- v^V/^ r wire) 
I 



Ji(-V^ly/JIawp) 



(B16) 



which describes a thin layer of excitation diminishing 
rapidly in the interior of the wire. 

We now want to determine the relationship between 



J 



the electric field and the magnetic induction at the wire's 
surface. The magnetic induction at the surface is a direct 
function of the driving current inside the wire. 



E z (p) = -dA z {p)/dt = icuA z (p) = - 



iluB„ 



surface 



I 



M- 



-ly/fUTWp), 



(B17) 



B =pH, 



(B18) 



B = V X A = e (( 



dA p {p) 
dz 

I 



o>-^>, 



(B19) 



we get 



B$_ 
P 



1 dA z (p) 
P dp 



^surface 



pJ\{— v— Ty 'pcrujr. 



-Ji(-V-1\/Jjiawp) 



wires 



(B20) 



which is proportional to the total driving current in the 
wire. Consequently 



(B21) 



E z {p) _ \/^l^pHjJ a (- -^ly/paujp) 
H 4> (P) \foJ\ (- V^y/pv^P) 



The two Bcssel functions grow at the same rate for large 
arguments and their ratio — > i. Consequently we have 
the wire surface result 



^( r wirc) = 
-^V( r wire) 



7^ 



(i-O 

V2 



x real parameters 



(B22) 



Thus the electric field along the wire has a constant phase 
shift relative to the driving current. This formula justifies 
the use of the complex constant times the measured high 
frequency resistance in the transmission line equations - 
again for constant a. 

For the case of spatially varying o we simplify mat- 
ters by assuming a semi-infinite slab of material with the 
vacuum-matter interface at x = 0, extending along y and 
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z to infinity. Assume the applied magnetic induction is 
parallel to the y axis. This implies the vector potential 
only has a z component. Because skin depths are much 
smaller than wire radii for most problems of interest - a 
ID approach is reasonable. 



V 2 A z (x) = 



d 2 A z 

dx 2 



-iui/j,a(x)A z (x) 



dA z {x) 

V X A — e y — \x=o— -"surface 



A z (x —> oo) = 



(B23) 



(B24) 



(B25) 



The simplest solution is to assume a Taylor's series 
expansion for all the terms in the equation. To simplify 
further we assume that the spatially varying conductivity 
is a small perturbation added to a constant conductivity. 



a(x) = (T + ec?i(x), £«1, 



A z (x) = A (x) + eA^x) + e 2 A 2 (x) + ■ 



d 2 A z (x) _ d 2 A (x) c d 2 A 1 (x) 



dx 2 



dx 2 



dx 2 



(B26) 



(B27) 



(B28) 



Applying these series to the diffusion equation we sepa- 
rate out the various contributions to the different powers 
of e 



d 2 A Q (x) 
dx 2 



= —iuii<JoAo(x) 



(B29) 



d 2 A 1 {x) 



-^-^ = -uJniaoA^x) + (ti(x)Ao(x)) (B30) 



The solution to the Aq equation with the correct 
boundary behavior at x — oo is 



A (x) = Cie -V-^^ox 

(-y/-lLU^(T x) n 

n=0 



(B31) 



An important thing to notice is that the Aq(x) solu- 
tion multiplies ai(x) in the equation determining A\(x). 
Since Aq (x) is an exponentially decreasing function away 
from the surface, the variations in o\{x) far away from 
the surface become irrelevant to A\{x). This explains 
why modeling the skin effect losses as a general function 
of frequency is necessary for accuracy. As the frequency 
changes, the skin depth changes - thereby averaging the 
wire's conductivity over a different spatial scale. 



To solve the A\ equation we expand as follows 

oo 

u x {x) = s £ j s n x n (B32) 

n=0 



Ay{x) 



n=0 



(B33) 



d 2 A 1 (x) ^ 

dx 2 = 2^n{n-l)a n x 



71=0 



= £(n + 2)(n + l) 



n=0 



(B34) 



CTl (x)A (x) 



oo n 

-EE 

n=0 j=0 



(-y/-lUJfJ,(To) n 1 

(n-j)\ 



(B35) 



Expanding the equations and separating by powers of x 
we find an infinite series of equations for the progressively 
higher order coefficients. The first two equations are 



= 



(B36) 



3 x 2a 3 + iuj^<T ai + ilo\ic\ (— y/— iujh<t ) x s 

+iwfj,c\ x si = 

(B37) 

Clearly ao and a\ are unconstrained by the recursion 
relationship. In the case where all the s n — then a n = 
as well - there being only a zero solution when there is 
no driving perturbation. Setting the first two coefficients 
to zero to match this limiting case, we see that the per- 
turbed solution only has finite terms of quadratic and 
higher power. This implies that at the wire's surface 
(x = 0) there is no change in either the vector potential 
or its curl from the perturbation in a. Thus the con- 
stant phase relationship between the electric field and 
the driving current in the wire remains as in the con- 
stant conductivity case. 
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